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Provided that matter-antimatter domains are present at the onset of big bang nucleosynthesis 
(BBN), the number of allowed additional relativistic species increases, compared to the standard 
scenario when matter-antimatter domains are absent. The extra relativistic species may take the 
form of massless fermions or even massless bosons, like relic gravitons. The number of additional 
degrees of freedom compatible with BBN depends, in this framework, upon the typical scale of the 
domains and the antimatter fraction. Since the presence of matter-antimatter domains allows a 
reduction of the neutron to proton ratio prior to the formation of 4 He, large amounts of radiation- 
like energy density are allowed. Our results are compared with other constraints on the number 
of supplementary relativistic degrees of freedom and with different nonstandard BBN scenarios 
where the reduction of the neutron to proton ratio occurs via a different physical mechanism. The 
implications of our considerations for the upper limits on stochastic gravitational waves backgrounds 
of cosmological origin are outlined. 
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I. INTRODUCTION 

The strongest constraint on additional energy density 
in the universe with a radiation-like equation of state is 
provided by big bang nucleosynthesis (BBN). The addi- 
tional energy density speeds up the expansion and cool- 
ing of the universe, and, consequently, the typical time 
scale of BBN is reduced in comparison with the stan- 
dard case. The additional radiation-like energy density 
may be attributed to some extra relativistic species whose 
statistics may be either bosonic or fermionic. Since the 
supplementary species may be fermionic, they have been 
customarily parametrized in terms of the effective num- 
ber of neutrino species 

N eG = 3 + AN cS , (1) 

where AN e g = corresponds to the standard case with 
no extra energy density. The standard BBN (SBBN) 
results are in agreement with the observed abundances 
for N c g = 2-4, giving thus an upper limit AN e g < 1 

Extra radiation-like energy can be also constrained 
by other physical considerations not directly related to 
BBN. However, the bounds seem to be much weaker than 
the ones provided by BBN. For instance, the growth of 
density perturbations in the early universe is affected by 
the radiation/matter energy density ratio, which in the 
standard case falls below unity some time before recom- 
bination. The data on CMB anisotropy @] has now im- 
proved to the point where it can be used to provide a 
constraint ]6|-|ll"| . Hannestad Q obtains an upper bound 
N cS < 17 (95% confidence), or AiV efi < 14. Using 



large scale structure (LSS) data in connection with the 
CMB data the same author also obtained a lower limit 
-^cff ^ 1.5/2.5 (depending on the data set used). 

The BBN and CMB limits are complementary. Both 
measure the increase in the expansion rate of the uni- 
verse, but through completely different physical phe- 
nomena. They also look at different epochs in the his- 
tory of the universe, BBN to T = Oil MeV), CMB to 
T = 0(1 cV). Thus attempts to relax one of these limits 
are unlikely to affect the other. 

The BBN limits can be relaxed in nonstandard BBN 
(NSBBN) scenarios (see the extensive reviews Jl2| and 
||, or Q). For example, the 4 He yield can be affected 
by changing the electron neutrino spectrum, like in the 
model with active-sterile neutrino mixing |15| , or by z/ e 
degeneracy, which changes the amount of thermal v e . 
This latter scenario is called degenerate BBN (DBBN) 

In the present paper yet a different scenario will be 
studied, namely BBN with matter-antimatter domains. 
In short the idea is the following. Consider the situation 
where matter-antimatter fluctuations |ll|[nj , are present 
at the onset of BBN. The fluctuations will be described 
in terms of their typical scale and the relative amount 
of antimatter. In previous studies (2(]]2l|] it has been es- 
tablished that this type of scenario is rather effective in 
reducing the neutron to proton ratio prior to the forma- 
tion of the 4 He abundance. The physical reason behind 
this statement is the different diffusion scale of protons 
and neutrons which can travel into antimatter domains 
annihilating there. Since the presence of extra relativis- 
tic species typically causes an increase in the neutron to 
proton ratio, it seems plausible that in the case of anti- 
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matter BBN (ABBN) more relativistic degrees of freedom 
can be accommodated than in the case of standard BBN 
scenarios. 

The additional relativistic species may be also bosonic, 
for instance gravitons. This possible interpretation of 
the extra energy density at the BBN epoch was orig- 
inally invoked by Schwartsman |E| since the produc- 
tion of gravitons is a generic phenomenon in Friedmann- 
Robertson- Walker (FRW) universes Relic gravita- 

tional waves decouple below the Planck scale and their 
energy density is constrained as part of the other ex- 
tra relativistic species. The BBN bound constrains the 
maximal fraction of critical energy density present to- 
day in relic gravitons of cosmological origin. The ex- 
perimental efforts are converging towards the possibility 
of direct detection of stochastically distributed gravita- 
tional waves . Various resonant mass detectors are 
now operating at a typical frequency of the order of kHz 
(Allegr o |26| , Auriga @, Explorer il, Nautilus [0 
Niobe 
600 



particles are photons, electrons, positrons, and N„ 
species of neutrinos, giving 



) and four Michelson interferometers (GEO- 
,TAMA-400 ||, LIGO |I| and VIRGO ||) 
will soon be operating in a wide frequency band from a 
few Hz up to 10 kHz. It is interesting, in the context of 
ABBN, to elaborate on the possibility that the additional 
relativistic species are of gravitational origin. 

The plan of our paper is the following. In Section II 
some basic considerations on BBN with extra relativistic 
species will be presented. In Section III the main features 
of ABBN will be summarized in light of the possible pres- 
ence of additional relativistic degrees of freedom in the 
scenario. In Section IV the results of the analysis of the 
parameter space of ABBN will be reported. In Section 
V the findings of ABBN will be compared with the case 
of DBBN, while in Section VI the constraints on the ad- 
ditional relativistic species obtained in the framework of 
ABBN will be studied together with other constraints 
coming from different physical considerations. Section 
VII deals mainly with the implications of our results for 
the case of relic gravitational waves. Finally Section VIII 
contains our concluding remarks. 
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The neutrinos have decoupled before electron-positron 
annihilation so that they do not contribute to the en- 
tropy released in the annihilation. While they arc rela- 
tivistic, the neutrinos still retain an equilibrium energy 
distribution, but after the annihilation their temperature 
is lower, T v = (4/1 1) 1 ^! 7 . Thus 

g* = 2 + 7 -N v {^\ = 2 + 0.454V, = 3.36 (4) 

after electron-positron annihilation. 

If we now assume that there are some additional rela- 
tivistic degrees of freedom, which also have decoupled 
by the time of electron-positron annihilation, or just 
some additional component p x to energy density with 
a radiation-like equation of state, p x — p x /3, its effect 
on the energy density and expansion rate of the universe 
is the same as that of having some (perhaps a fractional 
number of) additional neutrino species. Thus its contri- 
bution can be represented by replacing N v with N e g = 
N v + AN c g in the above. Before electron-positron anni- 
hilation we have p x — (7 /8)AN c gp~f and after electron- 
positron annihilation we have p x = 0.227A7V e ff/o 7 . 

The present ratio of the CMB (T = 2.725 K) en- 
ergy density to the critical density is f2 7 = p-y/p c = 
2.47 x 10~ 5 /i~ 2 . If the extra energy density component 
has stayed radiation-like until today, its ratio to the crit- 
ical density, Q x , is given by 



h 2 n r = K z — = 5.61 x 10~ 6 AV 



err 



(5) 



i 2 VL x = h 2 ^ 

Pc 

todayQ, where h is the Hubble constant in units of 
100 km/s/Mpc. 



II. BBN WITH EXTRA ENERGY DENSITY 

A. Preliminaries 

In the radiation-dominated era the energy density of 
the universe is 

This equation defines the effective number g r of relativis- 
tic degrees of freedom . An (ultra)relativistic fermion 
species with two internal degrees of freedom and in ther- 
mal equilibrium contributes 2-7/8 = 7/4 = 1.75 to g*. 
Before neutrino decoupling the contributing relativistic 



* To be more precise, the neutrinos have not completely 
decoupled by the onset of electron-positron annihilation, so 
that some entropy does leak to the neutrino component. This 
effect, plus finite-temperature QED corrections, can be rep- 
resented by using N e s = 3.04 after e ± -annihilation in the 
AiV c ff = case 

Mangano et al. Q have consid- 
ered the effect of extra relativistic species on this. Their nu- 
merical results (for AN c fc = 0-1) can be approximated by 
AW = 3.0395+(l-0.0014)1.0074AA 7 eff after e^annihilation, 
where the part 3.0395 - 0.0014 x 1.0074AAW is in neu- 
trinos and 1.0074(1 - 0.00025AA 7 e ff)AA r e fi is in the extra 
species [Bq|. This assumes that the extra energy density 
represented by AN c s has completely decoupled before e 
annihilation. Here we have defined AN c s as N c g — 3 be- 
fore e ± -annihilation, i.e., AN e g is a constant parameter, 
whereas N e g changes in the e annihilation. Then we have 
p x = 0.2288(1 - 0.00025A7V cf f) AiV cff p 7 after the annihila- 
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B. Standard BBN and extra relativistic species 

The BBN isotope most sensitive]] to the expansion rate 
is 4 He. The 4 He yield depends on the number of avail- 
able neutrons, or on the number ratio n/p of neutrons 
and protons. This, in turn, is determined by the com- 
petition between the expansion rate and the rates of the 
weak reactions which convert neutrons into protons be- 
fore the onset of the strong nuclear reactions which pro- 
duce the isotopes. Extra energy density speeds up the 
expansion rate so that more neutrons remain available 
for He production. 

In the usual terminology, standard BBN (SBBN) has 
AiVcff = 0, or N cB = N u = 3. Then SBBN is a one- 
parameter theory, with the present baryon-to-photon ra- 
tio 



r^lO 10 ^- 



(6) 



the only free parameter. It is related to fib, the fraction 
of critical density in baryonic matter today, by 



fl b h 2 = 3.65 x 10" 



(7) 



For the purpose of the present discussion we expand 
the concept of SBBN by allowing the presence of extra 
radiation-like energy density, so that SBBN has two pa- 
rameters, r\ and iV e fj. 

Since Y p , the primordial 4 He mass fraction, is an in- 
creasing function of both N e g and r\ in SBBN, the effect 
of a larger N e g can be compensated by a lower r\ to keep 
Y p in the observed range. 

The effect of the expansion rate on the other isotopes, 
D, 3 He, and 7 Li can be divided into two parts: 

1. The higher n/p ratio leads to an increase in the 
yield of the other isotopes too. 

2. The yields of these other isotopes are also sensitive 
to the time scale while these isotopes are produced 
by the strong nuclear reactions. The effect of less 
time can be compensated by increasing the reaction 
rates by increasing the baryon density. 

These effects work in the same direction for D and for 
3 He, but (for 7710 ^ 3) in the opposite direction for 7 Li, 
for which the second effect is stronger. Thus the net 



tionand tt x = (l-0.00025AiV eff )5.653x 10" 6 /i" 2 AiV off today. 
Since antimatter annihilation effects in our ABBN scenario 
take place in part during the same period as e ± -annihilation, 
the quantity AN c s constrained by our calculations is interme- 
diate between AN e a and N e s — 3 after e ± -annihilation. We 
ignore these small effects in the present work. 

^The magnitude of the effect on an isotope is considered 
in terms of how accurately its primordial abundance can be 
determined from observations. 



effect on the other isotopes is to shift the agreement with 
observations towards higher r\. This soon leads to conflict 
with Y p , constraining AN e s- 

Copi, Schramm, and Olive jij concluded that N e g < 4 
is a conservative upper limit. With various observational 
constraints on the primordial abundances, Lisi, Sarkar, 
and Villante obtained 95% confidence limits in the 
range AT c ff = 2-4. Using (probably prematurely) tight 
observational constraints, D/H = 3.4 ± 0.25 x 10~ 5 and 
Yp = 0.244± 0.002, and the prior AT eff > 3.0, Buries et al. 
|g| get a tight upper limit iV c ff < 3.20 (95% confidence). 
Cyburt, Fields, and Olive Q consider various fixed val- 
ues of rj (looking forward to the future precision determi- 
nations of r\ from the CMB anisotropy) and obtain the 
95% confidence limits N e s < 3.6(3.0) for 7710 = 5.8 and 
^cff < 4.1(3.9) for r/ w = 2.4 with (without) the prior 
N eS > 3.0. 

To summarize, the obtained SBBN upper limits are 



AN cS < 0.2... 1.0 



(8) 



depending on the observational constraints chosen. 
Translated into an upper limit to Cl x they become 



h 2 n T < 1.1... 5.6 X 10" 



(9) 



III. 



BBN WITH MATTER-ANTIMATTER 
DOMAINS 



BBN in the presence of antimatter regions has been 
discussed in popi] ]. In this scenario the baryon-to- 
photon ratio is not only inhomogeneous but also not 
positive definite. In spite of this the Universe is not 
matter-antimatter symmetric and the antimatter regions 
are small enough to be completely annihilated well be- 
fore recombination, so that they escape the CMB spectral 
distortion bound |39|. This type of inhomogcncities can 
come from different baryogenesis scenarios p9J . 

The time when most of the annihilation takes place 
is determined by the size of the antimatter regions, and 
this, in its turn determines the nature of ABBN. Two dif- 
ferent cases can be distinguished, depending on whether 
the annihilation takes place before or after 4 Hc is formed 
in nucleosynthesis. The relevant physics is completely 
different in the two cases. 

The case of interest here is the one with the smaller 
distance scale, where the antimatter region is annihi- 
lated before significant amounts of 4 He is made. This 
corresponds to antimatter regions of typical radius 
between 10 5 m and than 10 7 m. (We give all distances 
as comoving at T — 1 keV.) If the antimatter regions are 
even smaller they annihilate before neutrino decoupling 
and have no effect on BBN. The matter and antimat- 
ter are mixed by (anti)neutron diffusion. More neutrons 
than protons are annihilated and thus the n/p ratio is 
reduced compared to SBBN. Thus we get less 4 He. 
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This reduction in Y p due to neutron annihilation can 
be used to cancel the increase in Y p due to a large AN c g. 
Since both work by affecting the n/p ratio, the n/p part 
of the effect on the other isotopes is cancelled also, and 
one is left with the effect of the speed-up on the strong 
reactions. This shifts the agreement with observed abun- 
dances towards larger 77 for large AN c g. 

For calculations (see [^J for details), we assume an ide- 
alized geometry, where the antimatter regions are spher- 
ical with comoving radius ta and have homogeneous an- 
timatter density equal to the matter density outside the 
region. With this idealization, the ABBN scenario intro- 
duces two new parameters, ta and R, where R is the anti- 
matter/matter ratio in the universe. We assume R < 1 so 
that matter is left over after the antimatter regions have 
annihilated. Together ta and R determine the number 
density of the antimatter regions (or the typical distance 
between neighboring antimatter regions). 

In Fig. [l] we show the dependence of Y p on ta and R 
for 7710 = 6 and N e g — 3. We get the maximum effect 
for r ~ 10 6 -10 7 m, when most of the annihilation takes 
place just before 4 He is made. For fixed ta the effect 
is roughly proportional to R. For r < 10 7 m, the effect 
on the other isotopes is also via the n/p ratio and thus 
relatively small. 

For r > 10 7 m, some of the annihilation takes place 
after 4 He formation, leading to significant production of 
3 He and D. Because of 3 He overproduction, this region 
is not interesting for allowing larger N e g. 
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FIG. 1. The yield of 4 He as a function of the antimat- 
ter/matter ratio R and the radius ta of the antimatter re- 
gions, for 7710 = 6.0. The SBBN result, which is approached 
in the lower left corner, is Y p — 0.2484. 



IV. ABBN AND ADDITIONAL RELATIVISTIC 
SPECIES 

The observational constraints on the primordial abun- 
dances of the light elements are a matter of some contro- 



versy |40| . We follow Kneller et al. || in adopting the 
relatively generous constraints 

0.23 < Y p < 0.25 (10) 
2 x 10~ 5 < D/H < 5 x 10~ 5 (11) 
1 x 10~ 10 < 7 Li/H< 4 x 10~ 10 . (12) 

We search for the region in the parameter space which 
satisfies these constraints. In SBBN these constraints 
lead to the upper limit N e g < 3.4. For N e g = 3 the same 
constraints limit 77 to the range 4.2 < r/io < 5.9. 

For ta < 10 7 m, the effect of the antimatter regions 
on nucleosynthesis is mainly due to the reduction of the 
n/p ratio. As shown below, the dependence of the light 
element yields on ta and R can be combined into a depen- 
dence on a single parameter, the n/p ratio at the onset of 
nucleosynthesis. This redundancy of parameters allows 
us to cover the whole phenomenology by varying three of 
the four parameters only. For most of our computations 
we kept the antimatter radius fixed at ta = 10 6 ' 9 m, 
in order to get the maximal reduction in n/p for a given 
value of R. This left us with three parameters 77, N e g, and 
R. We calculated the ABBN yields of light elements in 
this 3-dimensional parameter space and compared them 
to the observational constraints. 

In Fig. ^|we show slices of this 3-dimensional parameter 
space for fixed values of 77. For a given 77, the upper limit 
on N e g comes from a combination of the upper limit on 
D/H and the lower limit on 4 He. 

In Fig. [3] we show slices along another dimension, for 
fixed values of R. The allowed N e g increases with in- 
creasing R. Simultaneously the allowed range in 77 shifts 
to higher values. 

Fig. |^ illustrates the redundancy of parameters R and 
ta- We show the yields of 4 Hc and deuterium as a func- 
tion of ta and R, for fixed 7710 = 6.0 and N g = 12. The 
value of 7V e ff was chosen close to its upper limit for this 
value of 7710 (see Fig. ||b) . Reducing ta moves the allowed 
region to higher values of R, but the width of the region 
remains nearly the same. Thus the effect of a smaller r a 
can be compensated by a larger R. The contours of Y p 
and D/H run almost parallel, so that we get essentially 
the same results for different combinations of ta and R, 
which lie along these contours. This shows that the light 
element yields depend on a single function of R and ta- 

In Fig. H we show the projection of the allowed re- 
gion in the (n,N e g,R) space onto the (77, 7V e ff)-plane. 
We see that very large amounts of extra radiation en- 
ergy can be accommodated in ABBN. A large N e g re- 
quires a somewhat larger 77. For, e.g, 7710 = 6.0 we obtain 
3 ^ N c g < 12.5. Even values iV c ff > 20 cannot be ruled 
out on the basis of nucleosynthesis. These large values 
require 7710 > 7. 

If we accept the possibility of ABBN, N e g is no longer 
constrained by BBN, but is limited by other observa- 
tional constraints. 
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FIG. 2. The light element yields in ABBN for a fixed 
r A = lO 6 ' 9 m and 7710 = 4.5 (a), 6.0 (b), and 9.0 (c), as a 
function of the two remaining parameters, R and N e g. We 
show contours for Y p (thick solid lines), log 10 D/H (thin solid 
lines), and log 10 7 Li/H (dashed lines). The shaded region is 
allowed by our adopted observational constraints. 
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FIG. 3. The light element yields in ABBN for r A = 10 6 9 m 
and antimatter/matter ratio R = 10 -1 ' 2 (a), 10 -1 ' 5 (b), and 
10~ 2 (c), as a function of 77 and N e g. The line types and 
shading are the same as in Fig. |ij 
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V. COMPARING ABBN TO DBBN 



A. Degenerate BBN 
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FIG. 4. The effect of the antimatter radius ta on the pro- 
duction of 4 He and D. We show contours of Y p (solid lines) and 
log 10 D/H (dashed lines) on the (rA, i?)-plane, for r/io — 6.0 
and iVefl = 12. The contours are almost parallel, showing 
that the yields depend on R and rA only through a single 
combination. 



It is interesting to compare these results to another NS- 
BBN scenario, degenerate BBN (DBBN) $l($\J§, which 
is also very effective in reducing the 4 He yield. In SBBN 
the neutrino chemical potentials are assumed to be small, 
so that there are almost equal numbers of neutrinos and 
antineutrinos. In DBBN this assumption is lifted. The 
degeneracy parameters 



S i rp 1 



(13) 



stay constant as the temperature falls, as long as the neu- 
trinos are relativistic. A nonzero increases the energy 
density in neutrinos of type i by 



30 
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AiV cff = ^ + -Jf. 



(14) 



Degenerate neutrinos are thus one possibility for a 
nonzero (positive) AN c g. 

If the degeneracy is in electron neutrinos, a much 
stronger effect is that it shifts the balance of the weak 
reactions between neutrons and protons. The equilib- 
rium ratio is 



Vl \ = e -Q/T-£° 

P/cq 



(15) 




FIG. 5. The combined allowed region in (r/,N c g). The 
shaded region indicates values of n, Ngg which are allowed 
for some value of the antimatter/matter ratio R. The dashed 
line is the CMB+SNIa constraint from fief. 



where Q = m n — m„. Thus a positive £ e reduces n/p and 
thus Y p , while a negative £ e increases n/p and Y p . This 
effect is much larger than the speed-up effect from the 
AiV e ff due to £ e - One can have both AN e g and £ e as free 
parameters in this scenario by allowing the other to 
have much larger (absolute) values than £ e (this possibil- 
ity might, however, be excluded by neutrino oscillations 
||41f ) , or by assuming some other extra radiation-like en- 
ergy density contribution. 

In DBBN one can compensate the increase in Y p due 
to a large AN e g by a reduction due to a positive £ e . Like 
in ABBN, one then compensates for the speed-up effect 
on the other isotopes by increasing rj. However this effect 
on these other isotopes is relatively small, and one can 
go to quite large AiV e ff with a relatively small shift in n. 

Considering BBN alone, the limit to AN e g is practi- 
cally removed. Combining DBBN limits with other cos- 
mological bounds on the density and expansion rate of 
the universe, Kang and Steigman obtained (in 1992) 
the limits —0.06 < £ e < 1.1, |^, r | < 6.9, which corre- 
spond to AJVeff < 71, and 2.8 < 7710 < 19. The other 
cosmological constraints have since gotten tighter (see 
Sec. M below). 
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B. ABBN versus DBBN 

The DBBN allowed region in r\ and AN e g (see, e.g., 
Fig. 2 of §, Fig. 5 of §], or Fig. 1 of @) is remarkably 
similar to that in ABBN (Fig. |). We can understand 
this as follows. 

In both cases the nonstandard feature affects nucle- 
osynthesis mainly by reducing the n/p ratio, with practi- 
cally arbitrarily large reductions possible. Thus the non- 
standard parameters, (r, R) in the case of ABBN and £ e 
in the case of DBBN, can be mapped onto a reduction 
factor in n/p. The rest of the game is then the same: 
A7V e ff can now be increased to bring the n/p ratio back 
up to give the observed value of Y p ; this has then a small 
effect on the other isotopes by shortening the time avail- 
able for nuclear reactions, which has to be compensated 
by an increase in r\. 



VI. DIFFERENT CONSTRAINTS ON 
ADDITIONAL RELATIVISTIC SPECIES 

A large radiation energy density delays (i.e., moves to 
lower temperature) the transition from radiation domi- 
nation to matter domination in the expansion law of the 
universe, and shortens the time scale of the radiation- 
dominated epoch. Since the radiation-dominated expan- 
sion law allows for logarithmic growth of density per- 
turbations only, the large-scale structure of the universe 
constrains the radiation component. The effect of the ra- 
diation component in the expansion law shows also as an 
integrated Sachs- Wolfe (ISW) effect in the CMB angular 
power spectrum. Hannestad 0] used this to obtain an 
upper limit N c g < 17 (95% confidence). 

In connection with a NSBBN scenario one can go fur- 
ther and combine the NSBBN constraints on (n,N e g) 
with the CMB constraints on (n,N c g). This has been 
done for DBBN in [@-|o|j4|| . These results will apply also 
to ABBN, since the allowed (n, N e g) region is essentially 
the same in both NSBBN scenarios. This NSBBN al- 
lowed region leans towards higher n for high iV G ff , whereas 
the CMB allowed region leans (slightly) towards lower rj 
for high N e g. The regions overlap for AN e g ~ few, but 
for the highest AJV cff allowed by CMB, the CMB seems 
to require a lower r\ than these NSBBN scenarios allow. 
Thus the combined constraint on AiV e ff may be lower 
than that from CMB alone. Hansen et al. jlO| use a dif- 
ferent set of CMB data, favoring somewhat lower n, than 
Kneller et al. || use, and get therefore tighter combined 
upper limits for AiV e ff. 

Since the CMB spectrum depends on a large number of 
other cosmological parameters, these results are sensitive 
to the assumed priors on these other parameters. Kneller 
et al. H find that values of AiV e ff larger than 10 are com- 
patible with combined CMB and DBBN constraints, but 
when they include a prior on the matter density parame- 
ter, < 0.4, and on the age of the universe, t > 11 Gyr, 



they get a tighter upper limit, A.N e g < 6. Hansen et al. 
|i~0f get AN c g < 5 with a prior on the age, t > 11 Gyr, 
and AAcff < 4 (both at 95% CL) with a prior on the 
parameters J7m and £Ia based on the supernova la data 
(see Fig. |). 

As emphasized by Hansen et al. [Q , the determination 
of n from the CMB data relies on the heights of both the 
second and the third acoustic peaks in the CMB angu- 
lar power spectrum, and, especially for the third peak, 
further data is needed to check on possible systematic 
effects. Therefore these combined CMB+NSBBN upper 
limits to AJVcff should be taken as preliminary. 

At present, we can conclude that values of AN e g as 
large as 4 are acceptable, and values somewhat larger, 
AN e g = O(10) may not yet have been ruled out. Thus, 
allowing for nonstandard BBN, like ABBN, raises the 
upper limit to the present fraction of extra radiation-like 
energy density from 



h 2 n T < 1.1... 5. 6 x 10" 



to 



h 2 n x < 2.3... 8 X 10~ 5 , 
an increase by an order of magnitude. 



(16) 



(17) 



VII. ABBN AND STOCHASTIC BACKGROUNDS 
OF RELIC GRAVITATIONAL WAVES 

The supplementary radiation allowed either by ABBN 
or by DBBN may be the consequence of a primordial 
background of gravitational radiation. The present crit- 
ical fraction of energy density stored in relic gravitons 



h 2 n 



GW 



( t „) = /j2PGW = 5 g x 1Q - 6 AiV c ff. 



(18) 



depends upon AiV c ff whose range of variation can be 
translated into constraints on the energy density of relic 
gravitational waves produced prior to BBN. The critical 
fraction of spectral energy density per logarithmic inter- 
val of frequency 



1 dp GW 
p r d In v 



(19) 




is often introduced for practical purposes, 
(obtained in the context of SBBN) and Eqs. 
(obtained in the context of ABBN) can then 
preted, in light of Eq. (|T8|), as an upper limit on the 
energy density presently stored in relic gravitons. Heed- 
ing experimental observations, the correlation of reso- 
nant mass detectors or of wide band interferometers can 
also provide upper bounds on the critical fraction of en- 
ergy density stored in relic gravitons. The sensitivity 
of operating resonant mass detectors is not yet able to 
probe possible signals compatible with Eqs. (p^|)-(^7|) 
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and ylq ) . Forthcoming generations of resonant mass and 
interferometric detectors are expected to improve their 
sensitivity in h 2 ilcw- The bounds coming from ABBN 
and from direct searches of stochastic gravitational wave 
backgrounds will represent two complementary sets of 
constraints on the same quantity. 



A. GW detectors and BBN bound 

Various resonant mass detectors are now operating 
p(i| po] | . In [[l3| , the first experiment of cross-correlation 
between two cryogenic detectors has been reported with 
the purpose of giving an upper limit on ft, 2 f2QW- The two 
detectors are Explorer |Q (operating in CERN, Geneva) 
and Nautilus [^9| (operating in Frascati, near Rome). 
Previous experiments giving upper limits on /i 2 r2cw used 
room temperature detectors. The Rome group obtained 
then an upper limit /i 2 J7gw < 60 at a frequency of 
roughly 905 Hz. The limit is a result of cross-correlation 
between the two detectors (located at a distance of ap- 
proximately 600 km) for an integration time of approxi- 
mately 12 hours. This limit is not competitive with the 
BBN bound (and also above the critical density bound 
implying that Oqw < !)• However, by increasing the 
correlation time from few hours to few months it is not 
unreasonable to go below one in /i 2 f2QW- 

Hollow spherical detectors have been recently inves- 
tigated as a possible tool for the analysis of the 
relic gravitational wave background. The sensitivity of 
two correlated spherical detectors could be 0(1O~ 6 ) in 
/i 2 r2cw for the frequency of resonance which lies between 
200 and 400 Hz. In this case the ABBN bound and the 
experiment will be certainly competitive. Dual spherical 
detectors may reach a sensitivity, in /i 2 f2QW, which 
is again O(10 -6 ) in the kHz region. 

Consider now the case of wide band interferometers 
§l]-|4| whose arms range from 400 m of TAMA up to 
the 3 km of VIRGO and to the 4 km of LIGO. The 
foreseen noise power spectra S n (v) are defined for v be- 
tween a few Hz and 10 kHz. While at small frequency the 
seismic noise dominates, at higher frequencies the main 
source of noise is provided by the shot noise. The signal- 
to-noise ratio (SNR) will be obtained by comparing the 
analytic form of the noise power spectra with the possible 
stochastic signal parametrized in terms of ncw(v). The 
single LIGO |34|] detector will be presumably sensitive to 
/i 2 J1gw ~ 10 . Two correlated interferometers lead to 
a SNR 



SNR 



3iJ ° 2 -FVTx 



2V2 IT 2 

7 2 (i/) tt GW (v) 



1/2 



(2), 



(20) 



(Hq is the present value of the Hubble parameter and 
F depends upon the geometry of the two detectors; in 



the case of the correlation between two interferometers 
F = 2/5). In Eq. ©, s£ k) (v) is the (one-sided) noise 
power spectrum of the k-th (k = 1, 2) detector, while 
7(2/) is the overlap reduction function [f46|-[49| which is 
determined by the relative locations and orientations of 
the two detectors. This function cuts off (effectively) 
the integrand at a frequency / ~ 1 /2d, where d is the 
separation between the two detectors. Eq. J20| ) assumes, 
as usually done, that the noises of the two detectors are 
stationary, Gaussian and not correlated. 

Consider now the simple case of a flat spectrum, 
namely Q,gw(v) — S} v°. From Eq. ( |2C| ) using the ap- 
propriate expressions of the noise power spectra and of 
the overlap reduction functions pertaining to the LIGO 
and VIRGO detectors j5(J we get that the sensitivity is 



h 2 n 



GW 



1.8 x 10" 



lyr 
T 



1/2 



SNR 



(21) 



for the correlation between two LIGO detectors [J49 50 
and 



h 2 n GW ~ 7.2 x 10" 



J SNR 2 (22) 



for the hypothetical correlation of two VIRGO detectors 
0. In Eqs. @ and @ T is the integration time of, 
for instance, one year. 

In the context of the standard BBN scenario we would 
be led to exclude h 2 ilow > 5.6 x 10~ 6 as a signal of 
primordial origin. This bound can be relaxed, according 
to Eqs. up to values which are O(10" 4 ), in the 

context of ABBN. If a signal with Ii 2 £Igw larger than the 
bound implied by standard BBN will ever be detected, it 
could still be of primordial nature. This observation may 
also be relevant for the bounds on stochastic GW back- 
grounds obtained by a single detector, where, h 2 flcw can 
be, at most O(10" 4 ). 

Microwave cavities have been originally proposed as 
GW detectors in the GHz-MHz region of the spectrum 
| pl| , p2| . Their application for high frequency gravita- 
tional wave backgrounds was theoretically suggested in 
p3[ . Improvements in the quality factors of the cavities 
(if compared with the prototypes of |5lJ| ) have been re- 
cently achieved |54j] and two experiments (in Italy |54| ] 
and in England R55| ) are now in progress. The sensitivity 
in /i 2 f2cw is still above one for the experiments reported 
so far. These devices could be important for the analysis 
of a frequency range where possible backgrounds gener- 
ated by collections of astrophysical sources are negligible. 



B. BBN bounds and theoretical models of stochastic 
GW backgrounds 

When the ordinary inflationary phase (of de Sitter or 
quasi-de Sitter type) is followed by a radiation epoch 
gravitational waves are produced and the theoretically 
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estimated f^cwM decreases as v~ 2 for 10 -18 kHz < 
v < 10~ 16 Hz. This branch of the spectrum corre- 
sponds to modes leaving the horizon during the inflation- 
ary phase and re-entering during the matter-dominated 
phase |5(]|37]]. Since f^cw^) is a decreasing function 
of the present frequency v the most relevant bounds 
(determining the normalization of the spectrum) will 
be the ones coming from small frequencies (i.e. large 
length scales) and, among them, the analysis of the first 
thirty multipoles of temperature anisotropics in the mi- 
crowave sky. The COBE normalization implies, in the 
low- frequency part of the inflationary spectrum, that 
h 2 fl GW (v) < 6.9 x 1CT 11 for v ~ 1(T 18 hKz. 

The modes re-entering during the radiation domi- 
nated e poch cor respond to frequencies 10 _16 Hz < v < 
10~ n y/Hi/Mp Hz, (where H x < KT 6 M P is the cur- 
vature scale at the end of inflation). The theoretical 
^Gw(ji) is either a flat (exact Harrison- Zeldovich spec- 
trum (58| ) or a slightly decreasing (in the quasi-de Sitter 
case) function of v, i.e. 0,Gw(p,to) cx v a , with a < 0. 
Imposing the COBE normalization as illustrated above, 
it turns out that ^gw(^) < 10" 15 for v > 10" 16 Hz. 
From Eqs. (EH])-© and |[|) it can be deduced that 
J]qw < 10" 1 ^vill always be compatible with the BBN 
constraints and out the foreseen sensitivity of wide band 
interferometers. 

In quintessential inflationary models ]59| ] , the inflation- 
ary phase is not immediately followed by a radiation- 
dominated phase, but a long stiff phase (dominated by 
the kinetic energy of the quintessence field) takes place 
prior to the dominance of radiation. As a consequence, 
^Gw(^) increases as v (up to logarithmic corrections) for 
10 _3 Hz < v < GHz Q. This branch of the spectrum 
corresponds to modes leaving the horizon during the in- 
flationary phase and re-entering during the stiff phase. 
Around the GHz, h 2 Q,Q^/{v) exhibits a spike. While the 
presence of the spike is a consequence of the fact that 
the inflationary epoch is followed by a stiff phase, the 
height of the spike is bounded by the BBN. In the case 
of standard BBN the height of the spike turns out to be, 
at most, 0.8 x 10~ 6 [|o|. Using the ABBN constraint the 
height of the spike becomes 1.7 x 10~ 5 . In pre-big bang 
models a similar discussion applies since the spectra of 
relic gravitons increase with frequency | |6l[ | and the most 
significant phenomenological bounds will be the one pro- 
vided by Eq. (|l8|). 

VIII. CONCLUSIONS 

Big bang nucleosynthesis provides the tightest con- 
straint on the radiation-like energy density in the uni- 
verse. However, some nonstandard BBN scenarios re- 
move this constraint. In particular, this is true for 
NSBBN scenarios which cause a large reduction of the 
neutron-to-proton ratio before 4 He is formed, such as 
ABBN and DBBN. 



All nonstandard features, whose effect on BBN is en- 
tirely due to a reduction of n/p before the nuclei are 
produced, have the same effect on BBN constraints on 
AJVeff, provided a sufficiently large n/p reduction is pos- 
sible. One can then consider just this n/p reduction as 
the nonstandard feature, and ignore the rest of the non- 
standard physics as far as BBN is concerned. 

The BBN upper limit to AN e g is raised, since this n/p 
reduction can be cancelled by a corresponding n/p in- 
crease due to the shortening of the timescale caused by 
a positive AN c g. This brings Y p back to the observed 
value. A small effect on the other isotopes remains since 
the timescale for the nuclear reactions is also shortened 
by AN c g. This can be compensated by an increase in n. 
Thus the larger AN c g values require somewhat larger val- 
ues of n. The required increase is about A7710 ~ A7V e ff/4. 

Allowing for nonstandard BBN removes the upper 
limit to additional energy density which comes from nu- 
cleosynthesis. Even when we include the other cosmo- 
logical constraints, the limit is raised from AN e g ~ 1 
at least to AN c g ~ 4, or from h 2 Q x = few x 10~ 6 to 
h 2 Q x = few x 10 -5 . Even values up to h 2 Q x = 8 x 10~ D 
may be acceptable. Since, after BBN, gravitational waves 
interact very weakly with the surrounding plasma, the 
upper limit of BBN sets, today, an upper limit on the 
critical fraction of energy density stored in gravitational 
waves of primordial origin. In the near future various 
detectors will be able to reach sensitivities, in /i 2 J7gWj 
which could be 0(1O -5 ) or even smaller. Hence, the 
problem will be to understand whether the obtained sig- 
nal is primordial or not. The bound on extra relativistic 
species obtained in the context of ABBN may be then 
relevant since ft, 2 f2cw cannot certainly be larger than 
10~ 4 but it can be larger than ~5x 10~ 6 , constraint of 
the standard BBN scenario. If a signal compatible with 
/! 2 Ugw > 5.6 x 10~ 6 will ever be found (either by two 
correlated resonant mass detectors or, most probably, by 
two correlated Michelson interferometers), it could still 
be of primordial origin provided /i 2 J7gw < O(10~ 4 ). 
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